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Abstract 

Let F be a nonlinear Frechet differentiable map in a real Hilbert space. Con- 
dition sufficient for existence of a solution to the equation F(u) = is given, and 
a method (dynamical systems method, DSM) to calculate the solution as the limit 
of the solution to a Cauchy problem is justified under suitable assumptions. 
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In this paper Theorem 1 gives a method for proving existence of the solution to nonlinear 
equation F(u) = and for computing it. Our method (dynamical systems method: DSM) 
consists of solving a suitable Cauchy problem which has a global solution u(t) such that 
y := u(oo) does exist and F(y) = 0. Theorem 1 generalizes Theorem 1 in [1] and its 
proof is based on the idea in [1]. The result of Example 1 in Section 3 below was obtained 



^ ■ in [5]. The result of Example 5 was essentially obtained in [8], where the existence of 



the solution was assumed, while in Theorem 1 below the existence of the solution is 
proved. In Example 6 some of the results from [4] are obtained. In Remark 2.4 in [8] 
some of the examples we discuss in Section 3 were mentioned. In [1] DSM was developed 
for solving ill-posed operator equations with not necessarily monotone operators and for 
constructing convergent iterative methods for their solution. 

Let F be a nonlinear Frechet differentiable map in a real Hilbert space. Consider the 
equation 

F(u) = 0. (1.1) 



*key words: inequalities, nonlinear equations, dynamical systems method, evolution equations 
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Let u) be a map, continuous in the norm of H and Lipschitz with respect to u in 
the ball B := {u : \\u — Uq\\ < R, u G H}. Weaker conditions, which guarantee local 
existence and uniqueness of the solution to (1.6) below, would suffice. Assume that: 

(F'(u)<f>(t,u),F(u))<- gi (t)\\F(u)\\ 2 VueB, (1.2) 

and 

\m,u)\\<g 2 (t)\\F(u)\\ VueB, (1.3) 
where gj,j = 1, 2, are positive functions on R + := [0, oo), g 2 is continuous, g± G Lj oc (R + ), 



and 



g x dt = +oo, (1.4) 

G(t) := g 2 (t) exp(- f 9l ds) G L 1 (R + ). (1.5) 
Jo 

Remark: Sometimes the assumption (1.3) can be used in the following modified 
form: 

\mt,u)\\<g 2 (t)\\F(u)\\ b VueB, (1.3') 

where b > is a constant. The statements and proofs of Theorems 1-3 in Sections 1 and 
2 can be easily adjusted to this assumption. 
Consider the following Cauchy problem: 

du 

u = $(t,u), u(0)=u , u:=—. (1.6) 

dt 

Finally, assume that 

\\F(u )\\ / G(t)dt<R. (1.7) 



The above assumptions on F, $, g\ and g 2 hold throughout and are not repeated in the 
statement of Theorem 1 below. By global solution to (1.6) we mean a solution defined 
for all t > 0. 

Theorem 1. Under the above assumptions problem (1.6) has a global solution u(t), 
there exists the strong limit y := lim t ^ 00 -u(t) = u(oo), F{y) = 0, u{t) G B for all t > 0, 
and the following inequalities hold: 



/oo 
G( 



x)dx, 



and 



\\F(u(t))\\<\\F(u )\\eM- [ 9i(x)dx). (1.9) 

Jo 

Theorem 1 generalizes Theorem 1 in [1]. 

In Section 2 proof of Theorem 1 is given, and two other theorems are proved, and in 
Section 3 examples of applications are presented. In Section 4 a linear equation and in 
Section 5 a nonlinear equation are discussed. 
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2 Proof of Theorem 1 and additional results. 



The assumptions about $ imply local existence and uniqueness of the solution u(t) to 
(1.6). To prove global existence of u, it is sufficient to prove a uniform with respect to 
t bound on ||tt(£)||. Indeed, if the maximal interval of the existence of u(t) is finite, say 
[0, T), and $(£,«) is locally Lipschitz with respect to u, then ||w(t)|| — > oo as t — > T. Let 
g(t) := \ \F(u(t))\\. Since if is real, one uses (1.6) and (1.2) to get gg = (F'(u)u, F) < 
~9i(t)g 2 , so 9 < ~9i(t)g, an d integrating one gets (1.9), because g(0) = \\F(uq)\\. Using 
(1.3), (1.6) and (1.9), one gets: 

pt px 

\\u(t) - u(s)\\ < g(0) G(x)dx, G(x):=g 2 (x)exp(- gi(z)dz). (2.1) 

J s JO 

Because G G L 1 (i? + ), it follows from (2.1) that the limit y := lim t _ >0O w(t) = u(oo) exists, 
and y G B by (1.7). From (1.9) and the continuity of F one gets F(y) = 0, so y solves 
(1.1). Taking t — > oo and setting s = t in (2.1) yields estimate (1.8). The inclusion 
u(t) G B for alH > follows from (2.1) and (1.7). Theorem 1 is proved. □ 
If condition (1.2) is replaced by 

(F'$,F)<- gi (t)\\F\\ a , 0<a<2, (2.2) 

then the proof of Theorem 1 yields g l ~ a g < —g\{t), so 

< git) < [g 2 - a (0) - (2 - a) / 9l (s)ds]^ . (2.3) 

Jo 

If (1.4) holds, then (2.3) implies g(t) = for all t >T, where T is defined by the equation: 

g 2 - a (0) - (2 - a) f gi (s)ds = 0. (2.4) 
Jo 

Thus ||F(«(i))|| = for t > T. So, by (1.3), $ = for t > T. Thus, by (1.6), 
u(t) = u(T) for t > T. Therefore y := u(T) solves equation (1.1), F(y) = 0, and 
\\u(T) - u(0)\\ < \\F(uo)\\ / T g 2 ds. If \\F(u )\ \ £ g 2 ds < R, then u(t) G B for all t > 0. 
We have proved: 

Theorem 2. If (1.2) is replaced by (2.2), (1.4) holds, and \\F(u )\\ J Q g 2 ds < R, 
whereT is defined by (2.4), then equation (1.1) has a solution in B = {u : \ \u—uq\ \ < R}, 
the solution u(t) to (1.6) exists for all t > 0, u(t) G B, u(t) := y fort > T, and F(y) = 0, 

yeB. 

If (2.2) holds with a > 2, and (1.4) holds, then a similar calculation yields: 

ft 

< g(t) < [^ (a ~ 2) (0) + (a -2) / gi (s)ds}^ := h(t) -> t -> 00, (2.5) 

Jo 

because of (1.4). Assume that 

poo 

/ g 2 (s)h(s)ds < R. (2.6) 
Jo 
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Then (1.3) and (1.6) yield ||u(f) -u(0)|| < R, \\u(t) -«(oo)|| < J t °° g 2 {s)h{s)ds -> as 
t — > oo. Therefore an analog of Theorem 1 is obtained: 

Theorem 3. If (2.2) holds with a > 2, (1-4) and (2.6) hold, then the solution u(t) 
to (1.6) exists for all t > 0, u(t) G B, there exists u{po) := y , and F(y) = 0, y G B. 

3 Applications. 

If gj = Cj, j = 1, 2, and Cj > are constants, then (1.4) and (1.5) hold, J °° Gdx = c 2 /ci, 
so condition (1.7) is: 

^||F(«o)||<i2, cj = Cj (R), j = 1,2. (3.1) 

Let us give examples of applications of Theorem 1 using its simplified version with gj = 
Cj > 0, j = 1,2. 

Example 1. Continuous Newton-type method. 

Let $ = — [F' (u)]^ 1 F (u) , and assume 

Ht^'H]" 1 !! < rri! = rm(R), Vm G B. (3.2) 

Assumption (3.2) holds in all the examples below. It implies "well-posedness" of equation 
(1.1). The case of ill-posed equation (1.1), when (3.2) is not valid, was studied in [2] for 
linear selfadjoint nonnegative-definite operator, in [3] for nonlinear monotone hemicon- 
tinuous operator, and in [1] for some other cases. Then c\ = 1, c 2 = mi, $ is locally 
Lipschitz if, for example, one assumes 

||F"(u)|| < M 2 , \/u G B, (3.3) 

where M 2 = M 2 (R) is a positive constant, and (3.1) is: 

mi (R)\\F(u )\\ < R. (3.4) 

In the examples below condition (3.3) is assumed and not repeated. As we have mentioned 
in Section 1, this condition can be weakened. 

Conclusion 1: By Theorem 1, inequality (3.4) implies existence of a solution y to 
(1.1) in B, global existence and uniqueness of the solution u(t) to (1-6), convergence of 
u(t) to y as t — > oo, and the error estimate (1.9). Condition (3.4) is always satisfied if 
equation (1.1) has a solution y and if Uq is chosen sufficiently close to y. 

Example 2. Continuous simple iterations method. 

Let $ = -F, and assume F'(u) > ci(R) > for all u G B. Then c 2 = 1, c x = ci(R), 
and (3.1) is: 

[dti^llFKlll <R- (3.5) 

If this inequality holds, then Conclusion 1 holds with (3.5) replacing (3.4). 
Example 3. Continuous gradient method. 



4 



Let $ = -[F']*F and assume (3.2). Then c 2 = M^R), because ||[F'(u)]*|| = 
\\F'(u)\\ < Mi(i2), and (F'$,F) = -| | [F'(u)]*F\ | 2 < -m^ 2 ||F|| 2 , so c x = m^ 2 , where 
mi is the constant from (3.2). Here we have used the estimates ||/|| = < 
II^IIIIA/H, \\Af\\ > p- 1 ]]- 1 !!/!!, with A := F'(u) and / = F(u). Estimate (3.1) is: 

M x m\ || F(u )|| < R- (3.6) 

If this inequality holds, then Conclusion 1 holds with (3.6) replacing (3.4). 
Example 4. Continuous Gauss-Newton method. 

Let $ = -([F'YF'y^F'YF, and assume (3.2). Then Cl = 1, c 2 = m 2 M 1; and (3.1) 

is: 

Mim?||F(u )|| < R. (3.7) 

If this inequality holds, then Conclusion 1 holds with (3.7) replacing (3.4). 
Example 5. Continuous modified Newton method. 

$ = — [F' (uq)}" 1 F (u) , and assume | |[F(m )] _1 | | < m . Then c 2 = m . To find c 1: let 
us note that: 

(F>$,F) = -IIFHH 2 - ((F'(u) - F'( Uo ))[F( Uo )]- 1 F 1 F) < 0.5||F( M )|| 2 , 
provided that 

\((F'(u) - F , (iio))[-F , (w )]" 1 i ? ,F)| < M 2 it>m ||F( M )|| 2 < 0.5||F( M )|| 2 , 

that is, M 2 i?mo < 0.5. If R = (2M 2 m ) _1 , then the last inequality becomes an equality. 
Choosing such R, one has c 2 = m , ci = 0.5, and (3.1) is: 2mo||i ? (uo)|| < (2M 2 mo) _1 , 
that is, 

4m 2 M 2 ||F(w )|| < 1. (3.8) 

If this inequality holds, then Conclusion 1 holds with (3.8) replacing (3.4). 
Example 6. Descent methods. 

Let $ = —-[jfj^h, where / = f{u{t)) is a differentiable functional / : H — > [0, oo), 

and h is an element of H. One has / = (/', u) = —f. Thus / = /oe - *, where / := f(u ). 
Assume ||$|| < c 2 |/| 6 , b > 0. Then ||w|| < c 2 \fo\ b e~ bt . Therefore tt(oo) does exist, 
/(tt(oo)) = 0, and ||w(oo) — u(t)\\ < ce~ bt , c = const > 0. 

If h = /', and / = ||F( M )|| 2 , then f(u) = 2[F']*(u)F(u), $ = -p^pf , and (1.6) is 

a descent method. For this $ one has c\ = |, and c 2 = where mi is defined in (3.2). 
Condition (3.1) in this example is condition (3.4). 
If (3.4) holds, then Conclusion 1 holds. 

We have obtained some results from [4]. Our approach is more general than the one 
in [4], since the choices of / and h do not allow one, e.g., to obtain $ used in Example 5. 



5 



4 Remark about linear equations. 



The following result was proved in [2]: if equation Ay = f in a Hilbert space has a 
solution y, and A > is a linear selfadjoint operator, then the global solution u to the 
regularized Cauchy problem 

u = -Au - a(t)u + f, u(0) = m , (4.1) 

has a limit lim^ 00 -u(t) := -u(oo), and A(u(oo)) = f. In [2] u G H is arbitrary, a > is 
a continuously differentiable, monotonically decaying to zero as t — > oo, function on P + , 
J °° ctrft = +oo, and a _2 d — > as t — > oo. 

Unfortunately the author of [7] was not aware of the paper [2]. Some of the proofs in 
[7] are close to these in [2] . The author thanks Dr. Ya. Alber for pointing out references 
[2] and [3]. 

If a > 0, a < 0, and a _2 |o;| < c, where c = const, then a~ l {t) — a _1 (0) < ct, so 
a(t) > [ct + a" 1 ^)]' 1 and consequently J °° adt = +oo (cf [7]). Therefore the condition 
/ °° adt = +oo in [2] can be dropped. 

In this Section we give a new derivation of the result in [2] under weaker assumptions 
about a, and show that the regularization in (4.1) is not necessary. 

First, let us prove that the regularization in (4.1) is not necessary: the result holds 
with a — 0. Below — > denotes strong convergence in H . 

The solution to (4.1) with a = is u it) = U(t)u + £U(t - s)fds, where U(t) : = 
exp(— tA). If E\ is the resolution of the identity of the selfadjoint operator A, then 
U(t)u = J °° e~ tx dE x u — > Pu as t — > oo, where P is the operator of the orthogonal 
projection on N, and iV is the null-space of A. Also J * U (t-s)fds = J °°(l — e~ tx )dE\y — > 
y — as £ — > oo, by the dominated convergence theorem. Thus, tt(oo) = y — Py + Pu 
and A(«(oo)) = /. □ 

Consider now the case < a — > as t — > oo. If /i(t) := exp(J * a(s)rfs), and u solves 
(4.1), then 

u(t) = h- 1 (t)U(t)u + h- 1 (t) exp(-tA) / e sX h(s)ds\dE x y. (4.2) 

Jo Jo 

Using L'Hospital's rule one checks that 

, A [*e sX h(s)ds , Ae u /i(i) 

lim Jo — -V = lim VTYTT^ A , ^ ^ = 1 VA > 0. 4.3 

t-oo e* A /i(t) t-oo \e tx h(t) + e* A /i(t)a(t) v y 

From (4.2), (4.3), and the dominated convergence theorem, one gets u(oo) = y — Py. 
The first term on the right-hand side of (4.2) tends to zero as t — > oo (even if Pu ^ 
0), if h(oo) = oo. To apply the dominated convergence theorem, one checks that 

AJ ° te ~ ( fe(tT Ma)da = XI ° e ~Kt) t ~ 8)d8 - 1 for a11 1 > and a11 A > °' where the inec i uai - 

ity < h(t — s) < h(t), valid for s > 0, was used. □ 
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Our derivation uses less restrictive assumptions on a than in [2] and [7]: we do 
not assume differentiability of a, and the property lim^oo a~ 2 a' = 0. The property 
J °° adt = +oo, which is equivalent to h(oo) = oo, was used above only to prove that 
lim^oo /i~ 1 (t)?7(t)-uo = 0. If / °° adt := q < oo, then tt(oo) = y — Py + e~ q Pu , and 
Au(oo) = /, so that the basic conclusion holds without the assumption h(oo) = oo. 

Finally, let us prove a typical for ill-posed problems result: the rate of convergence 
u(t) — > y can be as slow as one wishes, it is not uniform with respect to /. Assume a — 0, 
but the proof is essentially the same for 0<a^0ast^oo. Assume that A > is 
compact, and Aipj = Xjtpj, (ipj,ip m ) = 5j m . Then (4.2) with y — y m :— ip m and uq = 
yields u(t) = <p m (l — e~* Am ). Thus u(oo) = y, but for any fixed T > 0, however large, 
one can find m such that ||w(T) — y m \ \ > 0.5, that is, convergence is not uniform with 
respect to /. 

5 Remark about nonlinear equations. 

In this Section we give a short and simple proof of the basic result in [3] , and close a gap 
in the proof in [3], where it is not explained why one can apply the L'Hospital rule the 
second time. 

The assumptions in [3] are: the operator A is monotone (possibly nonlinear) , defined 
on all of H , hemicontinuous, problem (4-1) has a unique global solution, equation A(y) = 
f has a solution, a(t) > decays monotonically to zero, lim^oo aa~ 2 = 0, and a is 
convex. 

We refer below to these assumptions as A3). If A3) hold, the basic result, proved 
in [3], is the existence of -u(oo) := y, and the relation A(y) = f. In [3], p. 184, under 
the additional assumption (1.24) from [3], the global existence of the solution to (4.1) is 
proved. Actually, the assumption about global existence of the solution to (4.1) can be 
dropped altogether: in [6], p. 99, it is proved that A3) (and even weaker assumptions) 
imply that problem (4.1) has a unique global solution. 

Let us give a short proof of the basic result from [3]. It is well known that A3) imply 
that the problem A(v a ) + av a — f = 0, for any fixed number a > 0, has a unique solution, 
there exists \im a ^ v a := y, A(y) = f, and \\y\\ < \\z\\, for any z G {z : A{z) = /}. 
Thus, for any small 5 > 0, one can find as such that ||i> a — y\\ < 5/2 for all a > as, 
lim^oa^ = 0. Let w := u — v a , where u solves (4.1) and v a does not depend on t. 
Then w = — [A(u) — A{y a ) + a(t){u — v a ) + (a(t) — a)v a ]. Multiply this by w, use the 
monotonicity of A, and let ||w|| := g. Then gg < —a(t)g 2 + c\a(t) — a\g, c = \\y\\. 
Indeed, multiply A{v a ) + av a — A(y) = by v a — y and use monotonicity of A to get 
a(v a , v a — y) < 0. Thus | \v a \ \ < \\y\\, so c = ||y||. 

Since a(t) is convex, one has \a(t) — a\ < \a(t)\(t a — t), where t a > t is defined by the 
equation a = a(t a ), lim^o^a = oo. Thus, gg < —a(t)g 2 + c|o;(t)|(t Q — t)g, and, taking 
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u(0) = v a , one gets 

g(t a ) < ce-ti aa{x)dx [ " el° a(x)dx \a(s)\(t a - s)ds. (5.1) 
Jo 

We prove below that 

\ima(t)eti a{s)ds = 00. (5.2) 

This allows one to apply twice L'Hospital's rule to the right-hand side of (5.1), and to 
get: \im a _ g(t a ) = lim^-^ = °- Now > IK*<*) _ y|l < \\u{t a ) -v a \\ + \\v a -y\\, 

\\ v a~ y\\ < <^/2, and ||u(i a )— v a \\ < 5/2, for sufficiently large t a . Since 5 > is arbitrarily 
small, it follows that lim^oo | \u(t) — y\ \ — 0. 

Let us prove (5.2). From our assumptions about a, it follows that for all sufficiently 
large t, one has — aa~ 2 < c, where < c < 1, so a(t) > (ci + t^b, where b := c _1 > 1, 
ci > is a constant, and e# o(s)ds > (c x + t) b . Thus, (5.2) holds. □ 

If one assumes additionally that A is Frechet differentiate, then the proof is shorter. 
Namely, let h(t) := \\A(u(t)) + a(t)u(t) - f\\ := ||^||. Then hh = -((A'(u(t)) + 
ct(t))ip,ip) < —a(t)h 2 , because A' > due to the monotonicity of A. Thus h(t) < 0(t), 
where (f>(i) := h(0)e~fo ads . As we proved in Section 4, the assumptions on a(t) imply 
&{t) > {c-\t + C2) _1 , where C\ and C2 are positive constants, and C\ can be chosen so that 
< Ci < 1, due to the assumption lim^oo aa~ 2 = 0. Therefore J °° <j)(t)dt < 00. From 
(4.1) one gets: < (pit). Because J °° (f>(t)dt < 00, it follows that u(oo) := y exists, 
and ||m(oo) — u(t)\\ < <j)(s)ds. Finally, A(y) = f because h(oo) = = — /||. 

Any choice of a, for which 4>{t)dt < 00, is sufficient for the above argument. □ 
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